New Critical Point Induced by the Axial Anomaly in Dense QCD 
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We study the interplay between chiral and diquark condensates within the framework of the 
Ginzburg-Landau free energy, and classify possible phase structures of two and three-flavor massless 
QCD. The QCD axial anomaly acts as an external field applied to the chiral condensate in a 
color superconductor and leads to a crossover between the broken chiral symmetry and the color 
superconducting phase, and, in particular, to a new critical point in the QCD phase diagram. 
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Quantum chromodynamics at finite temperature, T, 
and chemical potential, (i, exhibits a rich phase struc- 
ture, from the hadronic Nambu-Goldstone (NG) phase 
at low T and ji, to the quark-gluon plasma (QGP) at 
high T, and color superconductivity (CSC) [lj at high 
/x. The phase transition from the NG phase to the QGP 
is being studied experimentally in ultrarelativistic heavy 
ion collisions at RHIC, and will be in the future at the 
LHC 0. The transition from the NG to the CSC phase 
could also be relevant in neutron and quark stars. 

The properties of the phases of dense matter in QCD 
depend on an important interplay between two compet- 
ing phenomena: quark- antiquark pairing, characterized 
by a chiral condensate (qq), and quark-quark pairing, 
characterized by a diquark condensate (qq) 0. We in- 
vestigate here how this interplay determines the transi- 
tion between the NG and the CSC phases, via a model- 
independent Ginzburg-Landau (GL) approach in terms 
of two order-parameter fields, the chiral condensate $ 
and diquark condensate d. This phenomenon is interest- 
ing not only in its own right, but also to other systems, 
e.g. the interplay between magnetically ordered phases 
and metallic superconductivity 

We consider two simple but non-trivial cases, 3-flavor 
quark matter with equal numbers of massless up (u), 
down (d), and strange (s) quarks, and 2-flavor quark 
matter with equal numbers of massless u and d quarks. 
Chiral and color symmetries rather stringently constrain 
the possible couplings between $ and d. As we show, 
the $-d coupling induced by the axial anomaly leads to a 
crossover between the NG phase and the CSC phase and 
also to a new critical point in the QCD phase diagram. 
The former may be relevant to the question of continuity 
of the hadronic and quark matter The case of two 
light quarks (u and d) and one medium-heavy (s) quark 
in /^-equilibrium with charge neutrality will be reported 
elsewhere. 

The GL free energy in three spatial dimensions in 
massless three- flavor QCD is £1 = £l x + fid + fl x d- The 
color-singlet chiral field $y ~ —(q^ R q\) is described by 



the standard free-energy 
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where "Tr" and "det" are taken over the flavor indices, 
i and j. Under SU(3) L x SU(3) R x U(1) B x U(l) A ro- 
tations, the chiral field transforms as $ — > c 2ioi a Vl$V r 
where the phase a A is associated with the U(\)a rota- 
tion. The first three terms on the right of Eq. are 
invariant under this rotation; the fourth term, caused by 
the axial anomaly, breaks the U(1)a symmetry down to 
%2N f = Zq. We assume that Co > 0, a necessary con- 
dition for the rf mass to obey mli > for positive $. 
We also assume that the chiral phase transition is driven 
by ao changing sign. The cubic determinant term in f2 x 
makes the chiral transition with three flavors first order. 

We focus on Lorentz scalar diquarks belonging to the 
fundamental representations in color and flavor space: 
((q L ) J b C(q L ) k c ) ~ e abc e ijk [dl} ai and ((q R Y b C(q R ) k c ) ~ 
£abc£ijk[d R ]ai where i,j,k (a,b,c) are the flavor (color) 
indices, and C is the charge conjugation matrix. Un- 
der SU(3) L x 5*7(3)^ x U(1) B x U(1) A x SU(3) C = Q, 

d L -> 



d L R transforms as d R 



,HcB-a A )v L d L V^; then (d R d{) 
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The most general form of the GL free energy of the d 
field, to C(d 4 ), is 0,110: 

fl d = a Ti[d L d\ + d R d R ] 

+ ([Tr(d L d{)} 2 + [Tr(d R d R )]- 



+ 2 (Tr[(d L dl) 2 ] + Tv[(d R d R ) 2 ]) (2) 
+ 3 Tr[(d R dl)(d L d R )} +/3 4 Tr(d L dl)Tiid R d R ). 

This free energy is invariant under Q. We assume that 
the normal-CSC transition is driven by ag changing sign. 
Since d^^R carries baryon number, det d^^ terms are not 
allowed, unlike for 

The interaction free energy of the chiral and diquark 
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fields is, to fourth order, 

n xd = 71 Tr[(d R d{)$ + (d L d R )&] 

+Ai Tr[(d L 4)$$t + (d R d R )&$] 

+A 2 Tr[d L 4 + d«4] ■ Tr[$t$] (3) 

+A 3 (dct$ • TrKdidJj)* -1 ] + ^-c-) • 

The triple boson coupling ~ 71, which breaks the E/(1).a 
symmetry down to Zq, originates from axial anomaly. 
The remaining terms are fully invariant under Q [10|. 
Equations Q-© constitute the most general form of the 
GL free energy under the conditions that the phase tran- 
sition is not strongly first order and that the condensed 
phases are spatially homogeneous. 

We assume, for three flavors in the chiral limit, a fla- 
vor symmetric chiral condensate, $ = diag(er, a, a), and 
a color-flavor-locked (CFL) diquark condensate, g?l = 
—d R = diag(d, d, d), in which all flavors contribute 
equally to the J p = + pairing. Then the GL free-energy, 
the sum of Eqs. Q-©, reduces to 

- 1 d 2 a + \d 2 a 2 . (4) 

Because the cubic d 2 a and a 3 terms in Eq. (@J both 
arise from the axial anomaly, 7 and c are related micro- 
scopically. Indeed, it is straightforward to show from the 
instanton-induced six-fermion interaction, &eti } j(q R q r L ), 
that 7 has the same sign and the same order of magni- 
tude as c 0,0. Positive 7 (attraction) favors coexis- 
tence, er 7^ 0, d 7^ 0. Also, the d 2 cr-term acts to break 
chiral symmetry explicitly, implying that er 7^ may be 
realized for all baryon densities, as discussed below. 

As we can show from microscopic calculations in weak 
coupling QCD and in the NJL model, A > and 
(3 > 0. A non- vanishing a plays the role of an effec- 
tive mass for the quark field, reducing the density of 
states at the Fermi surface, and the pairing energy 0, 
an effect represented \d 2 a 2 > 0. Furthermore, we find 
A//3 ~ ln(A/T ( j)/(/i/r £ j) 2 , which is rather small for rea- 
sonable values of [A, T d (the critical temperature of the 
color superconductivity without the a-d coupling), and 
A (~ /i for weak coupling QCD, and in the NJL model, 
~ the spatial momentum cutoff). We consider here the 
first order chiral transition driven by c with positive &; 
the case with negative b in three flavors does not change 
the results qualitatively fl~3 ) . On the basis of the above 
discussion, we focus on the case 7 > 0, A > with b > 0. 

In principle, the system can have four possible phases: 
normal (NOR) with a = d = 0, CSC (ct=0, d ^ 0), NG 
(a ^ 0, d= 0), and coexistence (COE) (a ^ 0, d ^ 0). We 
locate the phase boundaries and the order of the phase 
transitions by comparing the free energies, 
rt csc \0,d), ft( NG V,0), and ft( COE V,d). 




FIG. 1: Phase structure in the three-flavor system with 7 > 
and A = 0. The phase boundaries with a first (second) order 
transition is denoted by a double (single) line. 

In the COE phase, it is useful to analyze the free energy 
in terms of the single variable, cr(d), or d(a), obtained by 
solving the stationarity condition, dil^ COE \a, d)/dd = 0. 

In the absence of a-d coupling (7 = A = 0), the four 
phases defined above are separated by a = (a line of 
second order transitions) and by a = 2c 2 /(9&) (a line 
of first order phase transitions). But with an attractive 
coupling, 7 > with A = 0, the phase structure under- 
goes four major modifications, as shown in Fig^ where 
the first and second order phase boundaries are shown 
by double and single lines respectively: (i) The area of 
the COE phase grows, since the d 2 a term lowers the free 
energy if both a and d are nonvanishing. (ii) The first 
order line between the CSC-like and the NG-like COE, 
originally located at a = 2c 2 /(96), terminates at a critical 
point, A 0. This behavior is expected since d 2 a acts as 
an external field for er, washing out the first order phase 
transition for sufficiently large 7 or d. The NG-like COE 
phase in Fig.^has larger a than the CSC-like COE phase 
across the boundary BA. (iii) The second order bound- 
ary originally at a = splits in two, a line going to the 
right from the critical end point B and a line going to 
the left from the point C. Since a changes discontinu- 
ously across the first order boundary CB, the d 2 a term, 
which acts as a mass term for d 2 , leads to different criti- 
cal temperatures for diquark condensation to the NG-like 
COE and CSC-like COE. (iv) For 7 > (y/j3/b)c/3, a tri- 
critical point D appears on the phase boundary between 
NG and COE. Then the point C, otherwise a critical end 
point, becomes a triple point. 

Table I locates the characteristic points, A-D, in the 
(a, a) plane. To show explicitly how the critical point A 
appears, we derive an effective free-energy D,^p[a,d{a)], 
using the stationarity condition, d 2 = 2(ja — a/2)/ (3: 

cf 2 ci* c b 

n 3F [a, d(a)\ = - — + a* a + —a 2 - -a 3 + -a 4 , (5) 

with a* = arf//3 and a* = a — 2^ 2 //3. Equation (J5J 
is valid for a > 01/(27). We eliminate the a 3 term in 
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Eq. JSJ, by introducing r = a — cj (36). Then the system 
becomes equivalent to an Ising ferromagnet in an external 
magnetic field. The point A in FigQ] corresponds to the 
second order critical point of this magnetic system. 

The Xd 2 (T 2 term, with positive A, modifies the coef- 
ficients a*, c and b in Eq. JSJ as a* — ► a* — 2a\/(3, 
c -► c - 67A//3 and b -► 6 - 4A 2 //3. Therefore, for 
A//3 <C 1, as suggested microscopically, the A term does 
not qualitatively change the phase diagram in Fig^ 

Let us discuss the axial anomaly-driven crossover from 
the point of view of chiral symmetry. The CSC phase 
with a CFL structure {d L d R = —di&g(d 2 ,d 2 ,d 2 ) 7^ 0) 
breaks chiral symmetry but preserves the Z4 discrete sub- 
group of U(1)a- On the other hand, in the COE phase, 
d L d' R = —di&g(d 2 ,d 2 ,d 2 ) ^ and $ = diag(er, a, er) lead 
to chiral symmetry breaking, preserving only Zi- The 
symmetry breaking pattern is different in the two phases. 
However, the 71-term in Eq. J3J has Z§ symmetry which 
contains Z2 (but not Z4) as a subgroup. Therefore, once 
the axial anomaly is present, the NG-like and CSC-like 
COE phases cannot be distinguished by symmetry and 
can be continuously connected. The NG and COE phases 
differ in the realization of U{1)b symmetry, and therefore 
their boundary is not smoothed out. 

We turn now to the massless two-flavor system (with 
infinite s quark mass). In this case, all chiral and diquark 
condensates with an s quark are suppressed; we write 
$ = diag(cr, a, 0) and d L = -d R = diag(0, 0, d). The 
latter is the two flavor color superconductivity, 2SC state. 
Due to this color-flavor structure, the cubic terms in a 
and d are identically zero, and the model reduces to: 
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V) +AdV.(6) 



Since in two flavor QCD at finite T and fj,, a tricritical 
point (a = b = 0) may exist at which the second order 
transition for b > turns into a first order transition for 
b < |l5|. we introduce a a e term with positive coeffi- 
cient / as the minimal extension of the model to stabilize 
the system. 



TABLE I: For three flavors, with 7 > and A = 0: loca- 
tions of A, the critical point, B, the critical end point, C, 
the triple point, and D, the tri-critical point, in Fig. The 
point D appears for 7 > (^/ f3/b)c/3. For two flavors with 
b < and A > 0: locations of E and F in Fig. |3 Here 
c = V'c 2 + I667 2 //? and b' = b - 4A 2 //3. 
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FIG. 2: Phase strucuture in the two-flavor system for 6 > 
and A > 0. Left: Case with the tetra-critical point {\\fbp > 
A > 0). The second order line between NG and COE (CSC 
and COE) is characterized by a = 2aX/b (a = a/3/(2A)). 
Right: Case with the bi-critical point (A > ^^/bf3). The 
first order line between NG and CSC is characterized by 
a = av/ (3/b. 
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FIG. 3: Phase strucuture of the two-flavor system with 6 < 
and A > 0. The secon d order line between NG and COE is 
a = A(6 - sjb 2 - 4/o)//. The first order line between CSC 
and COE is at a = £[a - 3(6 - 4A 2 //3) 2 /(16/)]. The first 
order line between NG and CSC is at a 2 = (3[6bfa + (6 2 - 
4/a) 3 / 2 -& 3 ]/(6/ 2 ). 



We first consider b > in Eq. ©, neglecting the a 6 
term in finding the qualitative phase structure |l6(. For 
A = 0, the boundaries of the four phases are characterized 
by second order lines at a — and a = with a tetra- 
critical point at a = a = 0. With the repulsive d 2 a 2 term 
(A > 0), the area of the coexistence phase decreases, as 
shown on the left of Fig. [21 For A > -^\/b{3, the coex- 
istence region disappears, a first order interface between 
CSC and NG appears at a = a-J ' ft/b, and a = a = be- 
comes a bi-critical point, as shown in the right of FigEl 

Next we consider b < 0; here the a 6 term plays an 
essential role. For A = 0, the four phases are separated 
by a second order line at a — and a first order line 
at a = 36 2 /(!6/)- With a repulsive d 2 a 2 term (A > 0) 
the coexistence phase shrinks and gradually fades away as 
A — > 00. Moreover, a new first order line between NG and 
CSC appears and grows as A increases. This situation is 
shown in Fig. The locations of the critical end points 
E and F are given in Table I. Such a phase structure was 
previously noted in which used the random matrix 
model and in Q which used the NJL model; our model- 
independent analysis is consistent with these results. 

The mapping of the phase diagrams in the (a, a) plane 
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to the (T, /i) plane is a dynamical question which we 
cannot address within the GL theory. Nevertheless, in 
Fig. 01 we draw a speculative phase structure with two 
light quarks and a medium-heavy quark. There are two 
critical points in the figure: one near the vertical axis 
driven by finite fi (as seen in the two- flavor case), and a 
new critical point near the horizontal axis driven by the 
axial anomaly (as seen in the three-flavor case). With 
decreasing strange quark mass, m s , the higher critical 
point approaches the vertical axis, while as m s increases 
the lower critical point approaches horizontal axis. The 
chiral transition is a crossover in the direction of both 
high T and high fi. Whether this scenario is realized or 
not should be eventually checked by first principles QCD 
simulations. 




FIG. 4: Schematic phase structure with two light quarks (up 
and down) and a medium heavy quark (strange) . The arrows 
show how the critical points and the phase boundaries move as 
the strange-quark mass increases toward the two-flavor limit. 

In summary, we find that the QCD axial anomaly acts 
as an external field for the chiral condensate under the 
influence of the diquark condensate. The first order chi- 
ral transition is changed to a crossover for a large diquark 
condensate, and a new critical point driven by the axial 
anomaly emerges in the QCD phase diagram. Precise 
location this point is a future task for phenomenologi- 
cal models and lattice QCD simulations. Our schematic 
phase diagram would be made more realistic by includ- 
ing effects, such as finite quark masses, charge neutral- 
ity, P equilibrium, and thermal gluon fluctuations |H.ll8j. 
Open questions include whether the new critical point 
would survive in an inhomogcncous Fuldc-Fcrrcll-Larkin- 
Ovchinnikov (FFLO) state, and how the COE phase at 
low T and /i, Fig01 is affected by quark confinement. 
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